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Abstract. 

Halo coronal mass ejections (HCMEs) originating from regions close to the center of 
the Sun are likely to be responsible for severe geomagnetic storms. It is important to 
predict geo-effectiveness of HCMEs using observations when they are still near the Sun. 
Unfortunately, coronagraphic observations do not provide true speeds of CMEs due to the 
projection effects. In the present paper, we present a new technique allowing estimate the 
space speed and approximate source location using projected speeds measured at different 
position angles for a given HCME (velocity asymmetry). We apply this technique to 
HCMEs observed during 2001-2002 and find that the improved speeds are better correlated 
with the travel times of HCMEs to Earth and with the magnitudes ensuing geomagnetic 
storms. 
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1. Introduction 

Halo coronal mass ejections (HCMEs) originating from regions close to the 
central meridian of the Sun and directed toward Earth cause the most severe 
geomagnetic storms (Gopalswamy, Yashiro, and Akiyama, 2007 and refer- 
ences therein). Using coronagraphic observations it is possible to determine 
occurrence rate of CMEs, direction of propagation in the plane of the sky, 
angular width, and speed {e.g. Kahler, 1992; Webb et al, 2000; St. Cyr et 
al., 2000, Gopalswamy, Lara, and Yashiro, 2003; Gopalswamy, 2004; Yashiro 
et al, 2004). It is well known that geoffective CMEs originate mostly within 
a heliographic longitude ±30° (Gopalswamy et al, 2000, 2001; Webb et 
al, 2000; Wang, Sheeley, and Andrews, 2002; Zhang et al, 2003) and that 
the initial speed of the CMEs is correlated with the Dst index (Tsurutani 
and Gonzalez, 1998; Srivastava and Venkatakrishnan, 2002, Yurchyshyn, 
Wang, and Abramenko, 2004). These studies were based on the plane of 
the sky speeds of CMEs without considering the projection effects. Re- 
cently, there have been several attempts to correct for the projection effects 
and determine the real CME parameters (Zhao, Plunkett, and Liu, 2002 
(ZPL); Michalek, Gopalswamy, and Yashiro, 2003 (MGY); Xie, Ofman, and 
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Lawrence, 2004 (XOL); Michalek, 2006 (M)). These techniques, which need 
special measurements in coronagraphic fields of view, were based on the 
assumptions that CMEs have cone shapes and propagate with constant 
speeds. Based on the cone model, Michalek et al. (2006) demonstrated 
that fast CMEs originating on the western hemisphere close to the disk 
center are most likely to cause the most severe geomagnetic storms (Dst < 
-150 nT), although there have been some east-hemisphere CMEs capable 
of causing such kind of storms (Gopalswamy et al, 2005a; Dal Lago et al., 
2006) . Recently Moon et al. (2005) introduced a new direction parameter in 
correlation with gco-cffcctivcncss of CMEs. It is defined as the ratio of the 
shortest to the largest distance of the CME front from the disk center for 
a given LASCO image. They demonstrated that this parameter is strongly 
correlated with the geo-effectivcncss of CMEs. In this paper we present a new 
but similar technique to calculate the true speed of HCMEs, based on the 
dift'erences in the projected speeds measured around the occulting disk. Our 
direction parameters is determined from several (4-5) LASCO images and 
should be more accurate in describing the asymmetry of HCMEs. We use 
the space speeds for the prediction of travel times (TT) of HCMEs to Earth 
and to determine magnitudes {Dgr index) of geomagnetic storms. Results 
obtained with this simple technique are compared with those from other 
methods (cone models). The paper is organized as follows: in Section 2, the 
technique used to determine the improved speeds is presented. In Section 3, 
the data considered for this study is described. In Section 4, we implement 
the improved speeds for space weather forecasting. Finally, conclusions are 
presented in Section 5. 



2. The technique for the speed estimation 



The sky-plane measurement describing the properties of CMEs, especially 
HCMEs, are subject to the projection effects. However, true speeds in the 
direction of Earth are needed for estimating the travel times of CMEs to 
Earth (Gopaswamy et al., 2001). This issue is crucial for space weather 
forecasting. We propose to get realistic speeds by determining an asym- 
metry ratio (ASR) in the projected speeds measured at different position 
angles around the occulting disk. The ASR is defined as the ratio between 
the maximum and minimum speeds measured at different position angles 
around the occulting disk for a given halo CME. A two-step procedure is 
carried out to obtain the ASR. First, using height-time plots measurements 
at every 15° the projected speeds were determined at 24 position angles 
(see Figure 1 for an illustration of the method). Next, from these values, we 
calculate running averages from four successive speeds around the occulting 
disk. This allowed us to get 24 average speeds around the occulting disk for 
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each event. Thus the ASR is the ratio between the maximum (< Vmax >) 
and minimum (< Vmin >) average seeds determined in the way described 
above. There are two reasons to use the above procedure in determining 
ASRs. First, CMEs sometimes are irregular with narrow and fast structures 
which can over-estimate the measured projected speeds. Second, CMEs can 
be faint and single measurements arc prone to large errors. Our procedure 
minimizes these errors. It is important to note that the ASR has a real 
physical meaning. It is commonly assumed that CMEs have a cone shape 
which was described in details in a few papers (Michalek, Gopalswamy 
and Yashiro, 2003, Michalek, 2006, Xie, Ofman and Lawrence, 2004). Thus 
this parameter obtained directly from coronagraphic observations should 
better indicate CME propagation directions than could be deduced from 
the associated flare locations. Such locations of CMEs are very useful to 
eliminate the projection effects and to estimate the radial speeds and gco- 
effectiveness of CMEs. Moon et al. (2005) demonstrated that a direction 
parameter obtained from the asymmetry in the 3D shape of CMEs has a 
relatively good correlation with the geo-effectiveness of CMEs. Many halo 
CMEs have irregular shapes in the LASCO field of view, so the procedure 
proposed by Moon et al. (2005) is not accurate for such events. To minimize 
errors we propose to determine the directional parameter in the velocity 
space form several (4-5) LASCO images. Following their consideration, we 
defined the ASR using the maximum ratio between these projected speeds 
for a given event. We assume that the improved speeds of CMEs could be 
expressed by the formula: 

Vimp =< Vmax > + < Vmax > /ASR (1) 

where 

ASR =< Vmax > / < Vrmn > ■ (1) 

From these equations, we see that Vimp could be expressed as a simple sum 
of Vmax and Vmin- If a given CME originates at the limb then Vmin = 
and the radial speed is exactly equal to Vmax- If the location of a given 
CME moves from the limb to the center of the Sun then Vmax decreases 
and Vmin increases. At the disk center, when HCMEs are symmetric, Vmin 
= 2 < Vmax >. It is not a perfect approximation of the radial speed but the 
only possible because other method do not work for symmetric HCMEs. The 
projection effects are mostly determined by the source locations and widths 
of CMEs. Unfortunately, in comparison to the cone models, our technique 
allows us to estimate only the locations of CMEs. 
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Figure 1. An example of HCME (20 Jan. 2001) with dark lines along which the height-time 
measurements were performed. 

3. Data 

A list of HCMEs studied in this paper is displayed in Table 1. We considered 
only front-side full HCMEs during a period from the beginning of 2001 until 
the end of 2002. We selected this limited period of time to get a represen- 
tative sample of HCMEs to test our new tchnique. In the considered period 
of time, 70 front-side fuh HCMEs were found in SOHO/LASCO catalog. 
One of them was too faint to perform necessary measurements so 69 front- 
side HCMEs are displayed in the Table 1. Using data from World Data 
Center Kyoto (http://swdcdb.kugi.kyoto-u.ac.jp) geomagnetic disturbances 
caused by these events were investigated. In order to find a relationship 
between HCMEs and magnetic disturbances a two step procedure was used. 
First, we found all geomagnetic storms, in the considered period of time, 
with DsT index < — 30nT (this high limit was chosen following Michalek 
et al. (2006)). We assume that the associated magnetic disturbance should 
start no later than 120 hours after the first appearance of a given event 
in LAS CO field of view and no sooner than the necessary travel time of 
a given CME to Earth calculated from the measured maximal projected 
speed. We related a given disturbance with a given HCME if they were 
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within the specified time range. 20 events from our list were not geoefFective 
{DsT > — 30nT). These HCMEs were slow {V < 900kni/s) or originated 
closer to the solar limb. By examining solar wind plasma data (from Solar 
Wind Experiment, Ogilvie et al, 1995) and interplanetary magnetic field 
data (from Magnetic Field Investigation (Wind/MFI) instrument, Lepping 
et al, 1995), we identified interplanetary shocks driven by the respective 
interplanetary CMEs (ICMEs) . Measuring the time when a HCME first ap- 
pears in the LASCOs field of view and the arrival time of the corresponding 
shock at Earth the travel time (TT) can be determined {e.g. Manoharan et 
al., 2004; Gopalswamy et al, 2005c). 

The results of our considerations are presented in Table 1. The first col- 
umn gives the date of the first appearance in the LASCO field of view. Next 
three columns displays parameters obtained from the projected speeds (< 
Vmax >, < ymin. > and ASR). In column (5) linear speeds from SOHO/LASCO 
catalog are presented. Location of flares associated with CMEs are shown in 
column (6). In column (7) the minimum values of Dst indices for geomag- 
netic disturbances caused by HCMEs are displayed. Finally, in column (8) 
the travel time (TT) of shock to the Earth are presented. ASR varies from 
1.11 to 7.59. The lowest values are for events close to the disk center. The 
highest values are for limb HCMEs. Vmax is close to the catalog speed (V) 
for most events. The difference arises from the fact that Vmax is obtained as 
an average over 4 PAs, while V is form a single position angle. 



4. Space weather application 

For space weather application it is crucial to predict, with good accuracy, 
the travel times (TT) and strengths {Dst) of geomagnetic disturbances. In 
the next two subsections, we consider these issues using the improved speeds 
and the ASR parameter. 

4.1. Prediction of the onset of geomagnetic disturbances 

Figure 2 shows the scatter plot of the plane of the sky speeds (from SOHO/LASCO 
catalog) versus travel time. Correlation coefficients are 0.68 for the western 
and 0.50 for the eastern events. The standard error in determination of the 
TT is 16 hours. For comparison, we present a similar plot except that we use 
the the improved speeds (Figure 3, left panel). The figures clearly show that 
the improved speeds are better correlated with the TT, the correlations coef- 
ficients more significant: 0.72 for the western and 0.74 for the eastern events. 
The standard error in TT is 12 hours. The linear correlation coefficients 
are almost the same as those for the asymmetric cone model (Michalek, 
Gopalswamy and Yashiro, 2007). For (M) model the standard error was 
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49 Geoeffective halo CMEs, 2001-2002 
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Figure 2. Scatter plot of the plane of the sky speeds versus the travel times (TT). Dia- 
mond symbols represent events originating from the western hemisphere and cross symbols 
represent events originating from the eastern hemisphere. The dot-dashed line is a second 
degree polynomial fit to all the data points. 



smaller by ~2 hours. To illustrate that our considerations are consistent with 
previous results we compare them with the ESA model (the continuous line, 
Gopalswamy et al. 2005b). For these plots we used only the 49 geoeffective 
(DsT ^ — 30nT) events. For comparison, in Figure 3 (right panel) we have 
added the three events (2000/07/14, 2003/10/28, 2003/10/29) of historical 
importance, represented by the dark diamonds (Gopalswamy et al, 2005b). 
In our list there are CMEs originating from close to the disk center and all 
the way to the limb. For arrival at Earth we need to consider Earth-directed 
speeds rather than the space speeds. We obtained the Earth-directed speeds 
for the 49 geoeffective events and plotted them in Figure 3 (right panel). We 
see that points agree well with the ESA model. 

4.2. Magnitudes of geomagnetic disturbances 

Magnitudes of geomagnetic storms depend not only on the speeds of CMEs 
but also on the location of source region on the solar disk {e.g. Moon et 
al, 2005; Gopalswamy, Yashiro and Akiyama, 2007). Michalek et al. (2006) 
found a significant correlation between Dst and V *j where the parameter 
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Figure 3. The scatter plots of the improved (loft panel) and Earth directed speeds (right 
panel) versus the HCME travel time (TT). Diamond and cross symbols represent events 
originating from the western and eastern hemispheres, respectively. The dot-dashed line 
(left panel) is a second order polynomial fit to all the data points. The continuous line (right 
panel) is the ESA model representation. The three additional dark diamonds (only on the 
right panel) show the HCMEs (2000/07/14, 2003/10/28 and 2003/10/29) of historical 
importance. 



7 is the angular distance of a given CME from the plane of the sky. This 
parameter decides which part of a HCME hits Earth. Events with small 7 
strike Earth with their flanks while those with large 7 hit Earth with their 
central parts. Figure 4 shows the scatter plots of the plane of the sky speeds 
multiplied by 7 versus Dst index. The parameter 7 was determined from 
location of associated flares. There is a slight correlation between {V * 7) 
and Dst- Correlation coefficients are approximately 0.49 for the western 
and 0.30 for eastern events, respectivelly. Locations of the associated flares 
could be different from locations of the associated CMEs. For comparison, 
Figure 5 shows similar plot but for the improved speeds and the ASR pa- 
rameters. We propose to use the ASR instead of 7. This parameter includes 
source locations of CMEs (not flares) and should be more useful. There is 
only one difference to 7. Larger ASR means a location closer to the limb 
(7 1/ASR). Figure 5 shows the scatter plot of (< Vmax > + < Vmax > 
/ASR) /ASR versus Dst index. In other worlds, it is the scatter plot of 
improved speeds ((< Vmax > + < Vmax > /ASR) = Vimp) multiphed by 
source locations [l/ASR ~ 7) versus Dst index. We used similar formula in 
our previous paper (Michalek et al., 2006). There is a significant correlation 
between Vimp /ASR and Dst- Correlation coefficients are almost two times 
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69 Front halo CMEs, 2001-2002 
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Figure 4- The scatter plot of the plane of the sky speeds multiplied by 7 versus Dst 
index. Diamond symbols represent events originating from the western hemisphere and 
cross symbols represent events originating from the eastern hemisphere. The solid line is 
a linear fit to all the data points, the dot-dashed line is a linear fit to eastern events, and 
the dashed line is a linear fit western events. 



larger than those obtained from the projected speeds: -0.65 for the western 
and -0.61 for eastern events, respectively. It is clear that the our technic 
could be very useful for space weather applications. For these plots (Figure 4 
and Figure 5), we used all HCMEs form Table 1, even the non-geoeffective 
ones {Dst > —30nT). These events generate false alarms. Non-geoeffective 
HCMEs are slow {V <900km s^^) or have source region closer to the solar 
limb (see also Gopalswamy, Yashiro and Akiyama, 2007). The correlation 
coeffcient are slightly less significant in comparison to these received with 
the asymmetric cone model (Michalek, Gopalswamy and Yashiro, 2007). In 
the case of the cone model the correlation coefficient for all events was -0.72. 



5. Summary 

HCMEs are likely to cause severe geomagnetic storms. Speed is one of the 
most important parameter characterizing geo-effectiveness of CMEs. Un- 
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Figure 5. The scatter plot of (< Vmax+ < Vmax > / ASR) /ASR versus Dst index. 
Diamond symbols represent events originating from the western hemisphere and cross 
symbols represent events originating from the eastern hemisphere. The solid line is a 
linear fit to all the data points, the dot-dashed line is a linear fit to the eastern events, 
and the dashed line is a linear fit the western events. 

fortunately, coronagraphic observations, which provide CMEs speeds, are 
subject to the projection effects. In this paper, we presented a new way to 
estimate the radial speeds of CMEs. For this purpose we introduced the 
asymmetry ratio (ASR), the parameter obtained from the maximum and 
minimum speeds measured at different position angles for a given event. 
This allowed us to estimate realistic radial speeds. In Figure 2 and 3 we 
compared the travel time (TT) prediction using the projected and improved 
speeds. Respective correlation coefficients clearly show that the improved 
speed allows us to predict the TT with better accuracy. The standard error 
is about 2 hours smaller for the improved speeds than for the projected 
speeds. A similar situation appears when we compare correlation between 
speeds and source locations of CMEs and the strength of magnetic storms 
{Dst index). The correlation coefficients are about two times bigger for 
improved speeds in comparison with the projected speeds (Figure 4 and 
5). It is important to compare our results to previous considerations. Xie 
et al. (2006) calculated absolute differences between predicted (using the 
ESA model, Gopalswamy et al. 2005b) and observed shock travel times for 
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the cone models (XOL, MGY, ZPL). They found that the mean errors for 
those models were: 6.5, 12.8 and 9.2 hours, respectively. Recently, Michalek 
et al. (2007) calculated the same, but for the asymmetric cone model (er- 
ror=8.4 hours). In the present considerations, the mean difference between 
predicted (using polynomial fit from Figure 2) and observed shock travel 
times is 10 hours, similar to the error from complicated cone models. Many 
authors demonstrated that the initial speeds of CMEs are correlated with 
the DsT index but the correlation coefficient were not significant because 
they applied the plane of the sky speeds. Using the MGY model, Michalek 
et al. (2006) showed that the correlation between the space speed of HCMEs 
and DgT index could be much more significant (correlation coefficient was 
~0.60). Recently, Michalek, Gopalswamy and Yashiro (2007) showed that 
for the asymmetric cone model the correlation coefficient for the wester 
events could be very significant (0.85). In the present study wc found that 
correlation coefficient is 0.65 for the western events. Because in the present 
study we performed more accurate measurements of the projected speeds we 
got better results in comparison to the (MGY) model. Unfortunately, wc do 
not involve the width of CMEs so we can only estimate the radial speeds of 
CMEs and our results are not as accurate as those of the asymmetric cone 
model (M). There are two important advantages of this method. Firstly, 
this technique is very fast. To get the ASR we need only to determine 
the height-time plots from coronagraphic observations. We used the average 
maximum and minimum speeds obtained from four successive measurements 
running around the occulting disk. Secondly, we do not apply any model or 
assumptions. The method also has some weak points. Faint HCMEs can 
not be used for the study, because it is difficult to get the height-time 
plots around the occulting disk. Fortunately poor events arc generally not 
geoeffective so they are not of immediate concern. In contrast to the cone 
models we can not estimate the widths of HCMEs. Neglecting this parameter 
may lead to overestimation of the radial speeds. Although we obtained very 
accurate results using very simple considerations it is important to note 
that CMEs have very complicated 3D structure (Cremades and Bothmer, 
2004) and more factors may have to be included when determining their 
geo-effectiveness . 



Acknowledgements 



Work supported by NASA LWS T&T and SR&T grants. Work done by GM 
was supported in part by MMS' IF through the grant N203 023 31/3055. SY 
was supported in part by NASA (NNG05GR03G). 



paper21astv.tex.tex; 2/02/2008; 17:41; p. 10 



The space weather forecast 



11 



References 

Cremades, H., Bothmer, V.: 2004, Astron. Astrophys. 422, 307. 

Dal Lago, A., Gonzales, W.D., Balmaceda, L.A., Vieira, L.E., Echer, E., et al: 2006, J. 
Geophys. Res. Ill, A07S14. 

Gopalswamy, N.; 2004, ASSL series, cd. G. Poletto and S. Suess, KLUWER/Boston, 
p.201. 

Gopalswamy, N., Lara, A., Lepping, R.P., Kaiser, M.L., Berdichevsky, D., St. Cyr, O.C.: 

2000, Geophys. Res. Lett. 27, 145. 
Gopalswamy, N., Yasliiro, S., Kaiser, R.A., Howard, R.A, Bougeret, J.-L.: 2001, J. 

Geophys. Res. 106, 292907. 
Gopalswamy, N., Lara, A., Yashiro, S.; 2003, Astrophys. J. 598, L63. 
Gopalswamy, N., Barbieri, L., Lu, G., Plunkett, S.P., Skoug, R.M.: 2005c, Geophys. Res. 

Lett. 32, L03S01. 

Gopalswamy, N., Yashiro, S., Liu, Y., Michalck, G., Vourlidas, A., Kaiser, M.L., Howard, 

R.A.: 2005b, J. Geophys. Res. 110, A09S15. 
Gopalswamy, N., Yashiro, S., Michalek, G., Xie, H., Lepping, R.P., Howard, R.A.: 2005a, 

Geophys. Res. Lett. 32, L12S09. 
Gopalswamy, N., Yashiro, S., Akiyama, S.: 2007, J. Geophys. Res. 112, A06112. 
Kahler S.W.: 1992, Annu. Rev. Astron. Astrophys. 30, 113. 

Lepping, R.P., Acuna, M.H., Burlaga, L.F., Farrell, W.M., Slavin, J. A., et al.: 1995, Space 
Science Rev. 71, 207. 

Manoharan, P.K., Gopalswamy, N., Yashiro, S., Lara, A., Michalek, G., Howard, R.A.: 

2004 J. Geophys. Res. 109, A06109. 
Michalck, G.: 2006, Solar Phys. 237, 101. 

Michalek, G., Gopalswamy, N., Lara, A., Yashiro, S.: 2006, Space Weather J. 4, Cite ID 
S10003. 

Michalck, G., Gopalswamy, N., Yashiro, S.: 2003, Astrophys. J. 584, 472. 

Michalek, G., Gopalswamy, N., Yashiro, S.: 2007, Solar Phys. 246, 399. 

Moon, Y.-L., Cho, K.-S., Dryer, M., Kim, Y.-H., Bong, S., Chae, J., Park, Y.D.: 2005, 

Astrophys. J. 624, 414. 
Ogilvie, K.W., Chornay, D.J., Pritzenreiter, R.J., Hunsaker, F., Keller, J., et al.: 1995, 

Space Science Rev. 71, 55. 
Srivastava, N., Vcnkatakrishnan, P.: 2002, Geophys. Res. Lett. 29, 101029. 
St. Cyr, O.C., Howard, R.A., Sheely, N.R., Plunkett, S.P., Michels, D.J., et al.: 2000, J. 

Geophys. Res. 105, 18169. 
Tsurutani B.T., Gonzalez W.D.: 1998, Geophys. Monogr. 98; Washington, DC,AGU, 77. 
Wang, Y.M., Sheelcy, N.R., Andrews, M.D.: 2002, J. Geophys. Res. 107, 1340. 
Webb, D.F., Clivcr, R.W.,Crooker, N.U., Cyr, O.C., Thompson, B.J.: 2000, J. Geophys. 

Res. 105, 7491. 

Xie H., Ofman, L., Lawrence, G.: 2004, J. Geophys. Res. 109, A03109. 

Xie, H., Gopalswamy, N., Ofman, L., St. Cyr, O.C., Michalek, G., Lara, A., et al.: 2006 

Space. Weather. J. 4, S10002. 
Yashiro, S., Gopalswamy, N., Michalek, G., St. Cyr, O.C., Plunkett, S.P., Rich, N.B., 

Howard, R.A.: 2004, J. Geophys. Res. 109, A07106. 
Yurchyshyn, V., Wang, H., Abramenko, V.: 2004, Space Weather J. 2, Issue 2, CitelD 

S02001. 

Zhang, J., Dcrc, K.P., Howard, R.A., Bothmer, V.: 2003, Astrophys. J. 582, 520. 
Zhao, X-R, Plunkett, S.P., Liu, W.: 2002, J. Geophys. Res. 107, 1223. 



paper21astv.tex.tex; 2/02/2008; 17:41; p. 11 



MICHALEK ET AL. 

Table . 

List of halo CMEs with the determined asymmetry ratio ASR 
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